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Abstract 

The thermodynamic stability of odd-frequency pairing states is investigated 
within an Ehashberg-type framework. We find the rigorous result that in 
the weak coupling limit a continuous transition from the normal state to a 
spatially homogeneous odd-in-w superconducting state is forbidden, irrespec- 
tive of details of the pairing interaction and of the spin symmetry of the 
gap function. For isotropic systems, it is shown that the inclusion of strong 
coupling corrections does not invalidate this result. We discuss a few sce- 
narios that might escape these thermodynamic constraints and permit stable 
odd-frequency pairing states. 
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I. INTRODUCTION 



The puzzling physical properties of cuprate and heavy-fermion superconductors have mo- 
tivated theoretical investigations of new classes of superconductivity which are characterized 
by unconventional symmetries of the gap function A(k, cu^). A prominent example is the 
odd-frequency superconducting state where A changes sign under inversion of frequency, 
originally considered by Berizinskii in a model for triplet superfluidity in ^He [T]]. This idea 
was picked up by Balatsky and Abrahams who proposed a new class of odd-in-tu singlet su- 
perconductors which are odd under parity as well due to the Pauli principle |^. A possible 
microscopic realization of odd-frequency triplet pairing has been discussed by Coleman et 
al. in the context of heavy-fermion compounds 

The physical properties of these unconventional pairing states are rather unusual. Be- 
cause the equal-time gap vanishes, the order parameter of the condensed state is related 
to the expectation value of a composite operator 0,0. In case of a singlet pairing, the 
quasiparticle spectrum turns out to be gapless However, some concern has recently 

been raised with regard to the stability of odd-frequency pairing. As argued by Coleman et 
al., a uniform s-wave triplet odd-in-cj state should generically exhibit a negative Meissner 
effect and is therefore unstable with respect to a spatial modulation of the order parameter 
0. Additional support for a negative Meissner effect comes for calculations of the Meissner 
kernel using normal and anomalous Green's functions, whereas a positive Meissner effect is 
found within the composite-operator description [Q. Furthermore, uniform odd-frequency 
singlet pairing seems to contradict the requirements of causality and stability 

In this paper we intend to shed additional light on the question of stability of an odd- 
in-cij superconductor. This is achieved by examining the pairing induced change in the 
thermodynamic potential. The main result is that under quite general circumstances a 
continuous transition from a normal state into a spatially homogeneous condensed state is 
ruled out on thermodynamic grounds. This suggests that odd-frequency pairing cannot be 
realized in a uniform ground state. 
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The derivation of the above mentioned statement is presented in the next section. After 
introducing the formahsm in Section IIA, we first consider the weak couphng hmit, where 
the main ideas involved in the derivation can be presented in a very simple form (Section 
IIB). However, studies of the strong-coupling gap equation show that odd- frequency solu- 
tions often involve substantial renormalizations in the normal channel 0,0] • It is therefore 
important to include strong coupling corrections in the stability analysis. In Section IIC, we 
demonstrate that, at least for isotropic systems, strong coupling corrections are ineffective 
in stabilizing a uniform odd-frequency state. Implications of these results are discussed in 
Section III. 

II. STABILITY ANALYSIS 
A. Formalism 

We work within the framework of Fermi liquid theory. The following assumptions are 
made: (i) no spin-orbit coupling, (ii) pairing occurs in a quasiparticle band described 
by spin-independent energies Ck (measured from to the chemical potential), (iii) the low- 
temperature behavior can be described by a spin- and frequency-dependent effective inter- 
action V{aiki,(T2k2\cr3k3,a4k4^) where ai are spin indices and we have used a four-vector 
notation, k = (k, ti;„)- Fermi statistics implies that V is antisymmetric with respect to the 
first two and last two (cr, k) index pairs, respectively. The only further restrictions imposed 
on the functional form of V come from the conservation of total spin, energy and momentum 
iJ2i=i ki = 0), and from the hermiticity of the Hamiltonian which implies 

V*{aiki, a2k2\a3k3, aiki) = V{(yjii, cr3hW2k2, crih) (1) 

with k = (k, —iOn)- 

For a spatially homogeneous pairing state, the anomalous Green's functions are defined 

as 
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Jo 

F^^,{k) = fdre^^-^ < T,^lJ\^iT)^pl,iO) >, (2) 
Jo 

which are related by F^a'{k) = F*,^(k). Fermi statistics requires Fo-o-'(fc) = —F(ji^{—k). 
Using a matrix notation, the Gor'kov equations are given by 

[{lUn - ek)l - W{k)]G{k) = 1 - mF{k), 
[{lUn + ek)l + W{~k)]F{k) = -4>\k)G{k). (3) 

These equations contain the normal and anomalous parts of the electron self-energy defined 
by 

W^-Ti<72(^) = 4 2V{aik,a4p\a3p,a2k)G„^„^{p) 

P pcr[icr4 

(pa^a^Ak) = \Y^ V{crik,a2 - k\a3p, - p)F„^^Ap). (4) 

As a consequence of the antisymmetry of V, the relation (f)^^i{k) = —(f)^i^(—k) holds. 

The following stability analysis of the odd-frequency superconducting state makes use 
of a general expression for the change in the thermodynamic potential due to the two-body 
interaction V |^ 

Sn = U'^Y.tr{{^cOn-e^)G{k)-l},. (5) 

The trace is taken over spin indices, and the index g indicates that the quantities in the 
bracket correspond to a system with a scaled interaction gV . This formula is not restricted 
to weak coupling and is applicable even in the case of a critical lower coupling strength. We 
assume that the interaction potential V allows a certain type of solution of the gap equations 
(0) and (^, and examine its thermodynamic stability with respect to the normal state. It 
is now useful to distinguish between weak and strong coupling. 
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B. Weak coupling 



The weak coupling approximation is obtained by neglecting the self-energy contributions 
in the normal channel, W = 0, in which case 6Q equals the difference between the normal 
and superconducting state. Combining (|) and (^), we obtain 

'^^=":^r0?^"^'^^^)^^^^^^ ^^^^ 

= / Yg 5 + '^^^ + cP\k)m]-'cPHk)m}^ . (6b) 

Let us now consider a continuous phase transition. For temperatures close to Tc we can 
neglect the term proportional to 0^ in the denominator 

:c,- r tr{0t(k, -a;„)0(k,^„)}g 

In this form the sign of 5Vt can be determined by solely relying on symmetry properties of 
the gap function. There are two distinct cases: 

1. 0(k, — (^„) = (f){\^,Un), i.e. (j) is even in lo. Because 0^(k, ci;„)0(k, a;„) is a positive 
definite matrix (for all g)^ it follows that < and the paired state is stable with 
respect to the normal state. 

2. 0(k, —uOn) = — 0(k, cjn), i.e. (p is odd in u. This implies > showing that in this 
case the condensed state does not fulfill the thermodynamic stability criterion. 

Consequently, a second-order transition to a spatially homogeneous, odd-frequency super- 
conducting state is forbidden. Furthermore, in the weak coupling limit, this conclusion does 
not depend on the details of the pairing interaction and is valid for singlet and triplet pairing. 

C. Strong coupling 

We now extend the above given arguments to incorporate strong coupling corrections. In 
this case, no derivation with a generality similar to the weak coupling limit exists. However, 



we show that the same conclusions still hold for the class of isotropic models. This class is 
defined by 

(i) isotropic quasiparticle energies, ek = e|k|5 and particle-hole symmetry; 

(ii) a general spin-dependent interaction V(12|34) = r(1234) — r(1243) with 



r(1234) 



0"3 (T2Cr4 



4l-fc3,fc4-fc2 (8) 



where 1 = (kiai), etc. , and a*, i = x,y,z, denote Pauli matrices. F*" and are 
functions symmetric in frequency and rotationally invariant in momentum space; 

(iii) negligible |k|-dependence of F and the self-energy in the vicinity of the Fermi surface. 

The same type of interaction has been used in the work of Abrahams et al. |^. 

For singlet pairing the normal self-energy is diagonal in spin-space, W{k) = iuJnZ{k)l, 
and the gap function A is given by 0(A;) = ia'^Z{k)A{k). Under these conditions, is 
equivalent to 



Sn = - ['^ ljy^O)ly [ ^liZjk, ujn) - 1) + Z(k, a;n)A(k, a;n)A*(k, -a;„) 

Jo Q \ P n. J 4:71 ./,.,2 _L Aft \A*(t \ 



, (9) 
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where A/'(0) is the density-of-states per spin at the Fermi energy, and k is a vector on the 
Fermi surface. The renormalization factor takes the form 

^"/^-^ 471 J^^ + A(p,a;„)A*(p,-^^) 



Here, 7 = 2A/'(0)(F*-^ + 3F'^) is the effective interaction entering the normal self-energy. 

The possibility of a second order transition from the normal state into the superconduct- 
ing state is determined by the sign of the second order term in an expansion of in the 
gap function A. Using (|^) to replace Z, one finds 

"^^9 /■ G^f^k A(k,^„)A*(k, -a;„) 



g \ 471 \uJr\ 

7V^ ^ f dQ^^ f dQp Un . . i^m 



nm J 47r 

j" A(k,a;„)A*(k, A(p, a;^)A*(p, -o;^^) \ | ^^^^ 



Because 7 is symmetric with respect to an exchange (koj^) ^ (pcjm), the last term vanishes 
when summed over frequencies, and we are left with 

The expression for 6Q2 is very similar to (|^), and application of the same reasoning as in the 
previous subsection shows that strong coupling corrections cannot stabilize odd-in-a; pairing 
for the class of models considered. This result is in agreement with the conclusions obtained 
by Dolgov and Losyakov 0. 

The present derivation is not restricted to singlet pairing, but holds for unitary triplet 
pairing, too, where the normal self-energy retains the same spin-diagonal form, and (j){k) = 
iZ{k){d{k) ■ d)ay with d{k) x d*(A;) = 0. In this case, the product A{k)A*{k) is simply 
replaced by the scalar product d{k) ■ d*{k) in (p|)-(|l2l). No similar derivation seems to exist, 
however, for a nonunitary triplet pairing (d(/c) x d*{k) 7^ 0). 




III. DISCUSSION 

The result obtained in the previous section emphasizes some constraints imposed by 
thermodynamics which must be fulfilled by any realization of these unconventional super- 
conducting states. However, our derivation should not be mistaken as a proof of the nonexis- 
tence of odd-frequency pairing in general. As mentioned before, it requires the applicability 
of an Eliashb erg-type approach. Moreover, there are certain scenarios which are not covered 
by our analysis, and may allow stable odd-in-cu states. 

First, it is possible that in a system with large anisotropy, strong coupling corrections are 
more effective in producing a stable uniform odd-frequency superconducting ground state. 
Second, our analysis concentrated on second order transitions only, but does not exclude a 
first order transition. This possibility is easiest seen for singlet pairing in the weak coupling 
limit (i.e. (p{k) = zcr^A(A;)). In the odd-in-a; case, the term in the denominator of (|6bD 
is negative, leading to 
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1 f'dg^i |A(fc)p 
pJo g ^W + el-\A{k)\^ 

If the sign in the denominator is changed for a sufficiently large A;-region, the pairing state 
might be stabilized. 

A third scenario is based on the idea of a spatial inhomogeneity of the pairing amplitude. 
Here, the Cooper pairs acquire a finite center-of-mass momentum q, which describes the 
modulation of the order parameter on the lattice. This type of pairing has recently been 
discussed in more detail in the context of heavy fermion superconductors, and is thought 
to be driven by the negative Meissner stiffness of the homogeneous odd-in-a; state, which 
favors a coiling of the order parameter phase [0,0. 

It is easy to connect this picture to the thermodynamic stability consideration. Let (5f2(q) 
be the potential difference between the normal state and a superconducting state with a 
finite total momentum q of the electron pairs. The Meissner stiffness is proportional to the 
coefficients of the gradient terms in a Ginzburg-Landau expansion of the free energy, and is 
therefore related to the second derivative of 5i7(q) with respect to q. The positive Meissner 
stiffness found for a homogeneous phase of an even gap superconductor thus corresponds 
to a minimum of (5f2(q) at q = 0. The results of the previous section, in conjunction with 
a negative Meissner stiffness obtained for the uniform odd-in-cu state, suggests that in this 
case 5n(q) is instead maximal at q = 0. Consequently, the global minimum of 5i7(q) occurs 
at a finite momentum qmm, which defines an odd superconducting state with a positive 
Meissner stiffness. This state is stable with respect to the normal state if 6Q{cirnin) < 0, 
a necessary condition which may be violated even for interactions attractive in the odd- 
frequency pairing channel (i.e. which allows odd-in-iu; solutions of the gap equations). This 
contrasts the even-in-c<; case, where a homogeneous phase is always stable for an attractive 
pairing potential. 

One should keep in mind that the given discussion solely addresses the thermodynamic 
stability of odd-in-a; states with respect to the normal state. States which are stable accord- 
ing to this criterion do not necessarily correspond to the ground state, because competing 
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even-frequency states might possess even lower free energies. This question can only be 
answered when a specific interaction potential is given. 

In summary, we have shown that the thermodynamic stability criterion puts severe re- 
strictions on the occurrence of a continuous phase transition to a spatially homogeneous 
odd-in-w superconducting state. This result strongly suggests that any realistic model of 
this unconventional pairing requires a spatial modulation of the order parameter on an 
atomic scale from the very beginning. 
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